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The spherical model for spins describes ferromagnetic phase transitions well, but it fails at low 
temperatures. A quantum version of the spherical model is proposed. It does not induce qualitative 
changes near the phase transition. However, it produces a physical low temperature behavior. The 
entropy is non-negative. Model parameters can be adapted to the description of real quantum spins. 
Several applications are discussed. Zero-temperature quantum phase transitions are analyzed for a 
ferromagnet and a spin glass in a transversal field. Their crossover exponents are presented. 
0530, 0550, 7510J,7510N 

The spherical model was introduced by the late Marc Kac, see [1] for a discussion on its history and the route to 
its solution, presented by Berlin and Kac in 1952. [2] It was meant to be a simple soluble example for spin systems, 
from which the behavior near the phase transition can be studied in a simple manner. A review of this subject was 
given by Joyce. [3] 

The model also has a drawback: the behavior in the low temperature phase is pathological. The specific heat is 
just constant, while the entropy goes to minus infinity as S ~ logT. These aspects even occur for free spherical spins 
in a magnetic field. For these reasons the practical use of the spherical model has been limited. 
ON We have wondered whether the low T behavior can be healed by modifying the spherical model. We have considered 
the definition of spherical spins in a standard thermal field theory. We report our findings here. 
We consider spins coupled ferromagnetically to nearest-neighbors in the presence of a field 
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The Si are m-component variables. We shall mainly consider m = 3-component vector spins. The spins are subject 
to the spherical constraint 
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Later we will specify a. The classical partition sum reads 
^3 Z = f DSe' 1371 (3) 

a J 

S where DS = 6(^2^ =1 S? — Na) ]^[.{c/S;/(27r) m / 2 }. For free spins (J = 0) in a field one finds the free energy 
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The magnetization is M = H//i. Variation with respect to the "chemical potential" \i yields the equation of state 
M 2 = a — mT j fi. It can be solved explicitly 

M= (5) 

The susceptibility, x = & /mT, follows a perfect Curie law. In fixed field the magnetization goes to ^/a as T —> 0. 
However, the entropy goes to — oo as (m/2) logT for T —> 0. This is in conflict with the third law of thermodynamics, 
but occurs always for classical vector spins. The specific heat C = TdS/dT goes to a constant. 

We wish to improve this low temperature behavior. If we return to eq. (1) at H = and denote the eigenvalues of 
the coupling matrix by A, we obtain Z ~ max^ exp(Am/3 / u/2) Yl\W(M~^)]~ m ^ 2 ■ ^ n analogy with harmonic oscillators, 
we expect that the term f3(ji — A) arises from a quantum expression sinha/3( / u — A). Such a factor can indeed be 
derived when no square root occurs. We therefore propose two ways to quantize the spherical model. If spins are real 
valued, the quantization must be defined for Z 2 . This is reminiscent of the solution of the 2-d Ising model, where Z 
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can be evaluated with Pfaffians, while Z 2 can be expressed in terms of Grassmann integrals. If the spins are complex, 
the quantization can be defined for Z itself. 

For real valued spins, one doubles its components. One thus considers spins S?(r) (a = 1,2) defined on an 
"imaginary time" interval < r < [3, subject to periodic boundary conditions S°(/3) = S"(0). In units where Ti = 1 
the square of the partition sum reads 
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Z 2 = /2>Sexp / dr\—J2 E SUr)^^S»(T)-J2mSUr)})\ (6) 
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where the r-integral is discretized in M = /3/dr steps and dS(r) = S(t + dr) — S(t). The two spins systems 
are coupled by the Pauli-matrix in the derivative term, viz. erf, 2 - 1 = — = Cri = ^22 = 0- The integration 
measure is a normalization constant times the repetition of the sperical measure at all imaginary time points, VS = 
C Nm Yla,r DS a (r) with C = (l/2a) EL^o^M/ 2 ")- 

The quantum term in (6) involves a first order derivative in time. Our motivation hereto was that rf-dimensional 
quantum spins systems, such as the Ising chain in a transverse field, involve such a term due to the Trotter-Suzuki 
transformation into a classical d + 1-dimensional problem. In a recent work on the spherical model other systems 
were considered where a second order derivative occurs. [4] For magnetic systems this does not lead to acceptable low 
temperature behavior. 

In order to fix some universal properties of the system, we consider again free spins in a field. At each r we 
introduce a Lagrange multiplier h(t) for the spherical constraint. The mean field value occurs for constant pt. If the 
thermodynamic limit N —> 00 is taken before the limit M —> 00, the fluctuations do not yield extensive contributions. 
We thus find simply, instead of eq. (4), 

F = -^- + — log2sinha^-— (7) 



Variation with respect to jj, now yields 
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M 2 = — — = a — ma coth a/3a (8) 
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For large T, [i grows linearly, /j, = (y/a)T with coth y = a/ma. The susceptibility follows the Curie-law x = V/" = 
a/yT. The entropy reads 

Tfl 

S = — (a/3/j, coth a/3/j, — log 2 sinh a/3/j,) (9) 

Whereas in the classical model there was no unique way to define its adjustable constant, we have chosen it here such 
that, in a finite field, S vanishes at T = 0. In the classical limit a ~ Ti 2 —> one might be tempted to neglect the 
term — (m/2)loga. It is exactly this procedure, however, that produces a negative classical entropy. 
In zero field, or at infinite T, the entropy becomes 

Soo = ^ log 2 sinh j/) (10) 

2 \ma / 

The zero-point magnetization Mo = \/a — ma exhibits a quantum reduction from the classical value yfa ' . Eqs. (7) 
and (9) exhibit an energy gap AE = 2a / uo = 2aH/Mo- This is indeed expected for free quantum spins. 

We now wish to see in how far our spherical model can be used to describe realistic quantum spins. In our definition 
of the thermal field theory we had two parameters: the length a and a ~ k 2 = 1. We can fix these parameters by 
comparing to two expressions from the high and/or low temperature thermodynamics of free spins. Let us therefore 
consider quantum spins of order S where S = 1/2, 1, 3/2, • • • is half-integer or integer. Our field H is the product 
of the external field Ho, the (/-factor and the Bohr magneton hb, viz. H = g^BHo- For free quantum spins the 
susceptibility follows a Curie law 

g 2 l, 2 B S(S+l) 

Xo = 9VbX = ^ (11) 

This law is reproduced by our model with m = 3 if tanh3a/5'(5' + 1) = 3a/c We can find a second constraint by 
adjusting the zero point magnetization to Mo = S. This implies a = S 2 + 3a. Combining these relations one has to 
solve 
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tanhy = + ^ (12) 

y S+(S+l)y K ' 

This equation has a solution for any positive 5. We find the values, presented in the table. For large 5, a is close to 
5(5 + 1) + 1/3, while a « 5/3+1/9. The energy gap, AE/H « 2/3 + 2/95, is below the value AE/H = 1 of spin-5 
quantum spins, but it is bounded by 2/3. 

In Figure 1 we present the specific heat of a free 5=1 spin with the values of a and a taken from the Table. It is 
compared to the "classical" result (same a; a — ► 0), which goes to 3/2 at T = 0. We have also presented the specific 
heat of a 5 = 1 quantum spin, viz. C = (/3H/2 sinh(/3Jf/2)) 2 — (3/3H/2 sinh(3/3Jf/2)) 2 . This expression vanishes 
quicker at small T, since it has a larger gap than our result for the free spherical spin. 

Having fixed the system parameters we can now consider coupled spins. The first case is the mean field ferromagnet, 
eq. (1) with the sum involving all pairs. We denote Jo = N J . The free energy reads 

F = + ^logSsmh^ + ^1 _ (JoM + Hl 

2 2 8 1 ^ 2 2/j K ' 

When M ^ in zero field, one has fj, = J and 

M 2 = a - am coth a(3J (14) 

At the phase transition M grows as (T c — T) 13 with critical exponent j3 = 1/2, as in the classical description. The 
entropy in the ordered phase 

5 = y [a/3J (coth a/3J - 1) - log (l - e - 2al3J °)] (15) 

vanishes exponentially at low T. A similar behavior occurs for Ising spins with mean field couplings, since the latter 
just act as an external field, of adjusted strength. 

It is known that the spherical spins behave differently in case of short range ferromagnetic interactions. In d- 
dimensions one finds 

aa mT f d d k , . , ^ H 2 , , 

F = -- + — J w? log 2 smh a(][ " - m] - W^) (16) 

where Jo = J(k = 0). Our quantum formulation brings no qualitative changes near the phase transition. One still 
has [3 = 1/2, r) = 0, while also the other exponents are unchanged. However, near the transition the low temperature 
branch of the specific heat is linear rather than constant. At low T the excitations are spin waves. Their degeneracy 
is m, and not m — 1, because, apart from the direction, also the length of the spins can fluctuate. The entropy and 
the specific heat decay as mT d l 2 at small T. In Figure 2 we present a plot of the specific heat for 5=1 spherical 
"Heisenberg" spins on a simple cubic lattice, with a and a taken from the table. It is compared to the classical result 
(a = 0). 

Another case of interest is when only the z-components of the spins have an exchange coupling. The field may have 
a component in the transverse direction; a well known example is the Ising chain in a transverse magnetic field. For 
m = 3-component spherical spins one finds the free energy 

„ utr T f d d k , . , 

F = -^- + -J ^-plog2sinha^- J(k)} 

+ Tlog2smh a/3 ,-^-^^ (17) 

In a perpendicular field it leads to the equation of state 

/d d k H 2 
^ ^ d coth a/3\p — J(k)] — 2a coth a(3/jL j- (18) 

For d > 2 a phase transition occurs where M z ^ 0. The ordering is suppressed when Jf j_ is large enough. At T = 
there occurs a quantum phase transition at critical field Jfj_ = Jo 5. Generally, there is an energy gap, 2aH z /M z ; it 
only vanishes in the ordered phase at zero field. There is a scaling region when A = Jo5 — Jfj_, H z , T, and f3H z /M z are 
small. Here the singular part of the free energy and the order parameter assume the scaling form F s i ng = A 2 $i(a;, y), 
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M z = A 1 / 2 $ 2 (*,J/) with x = H z /\A\' t ' H , y = T/\A\' t ' T with crossover exponents <f> H = 2/3, <f> T = 2/d. At non-zero 
H z and very small T the thermal excitations experience the gap. 

A further subject of interest is the limit of strong anisotropy for Heisenberg spins. For considering the uniaxial 
limit where the spins point in the z-direction (Ising-like) we add an anisotropy term 

Hani = D y ] S ix + Sfy (19) 

i 

to the Hamiltonian. For large D the free energy reduces to eq. (16) with m = 1, provided we subtract the 'zero point 
energy' 2aD and make the shift a —> a = a — 2a. Under these modifications we can fully neglect the orthogonal 
degrees of freedom. The zero point magnetization remains equal to Mo = S. When first D — ► oo, the T = oo entropy 
reduces to 

V 1 

Soo = — 5- - -log2sinh?/ (20) 
la I 

This is below the isotropic value, as it should. 

In systems with XY-spins, the z-component of the spin is suppressed. Adding a term D^^S^ — a) we find back 
eq. (16) with m = 2, provided we take a —> a = a — a. 

Next we investigate quantum effects in the mean field spherical spin glass. This system is described by eq. (1) 
where now for each pair the coupling J 8 j is an independent Gaussian random variable with average zero and variance 
J 2 /N. The model was introduced by Kosterlitz, Thouless and Jones [6], and extended to the situation with short 
range ferromagnetic interactions by the present author [7]. In zero field the quantum free energy is given by eq. 
(16), provided we replace J d d k/(2iT) d by J dJ^AJ 2 — J 2 /2irJ 2 , involving the semi-circular eigenvalue density of the 
random coupling matrix. There is a phase transition where \i sticks at 2J 2 . Here the Edwards- Anderson parameter 
q becomes non-zero. At low T the entropy and specific heat now vanish as T 3 I 2 . In this mean field model there is no 
gap. 

For a spherical spin glass with couplings only in the z-direction, we can again look at the effect of a perpendicular 
field, while H z = 0. We find a similar analog of eq. (17), while in eq. (18) q should replace M 2 . At T = there is 
a quantum phase transition for if j_ = JoS, where the spin glass order is suppresed. In the disordered phase A < 0, 
there is again a gap. In the regions where the gap vanishes or can be neglected, the singular part of the free energy 
and the order parameter have the scaling form F, ing = A 2 $!(T/| A| 2 / 3 ), q = A$ 2 (T/| A| 2 / 3 ). 

For complex valued spherical spins we consider the thermal partition sum 

Z = |2)Sexp^d r |zl^ S *( r )i :S . (r) _ W(S . (r)) | (21) 

with VS = C N Y[ T DS(t) and DS = 6(2Na — S* Si)Y\i(dS* dSi/2ir) being the spherical measure for complex 
spins. The Hamiltonian now reads 



H 



■\jY, S *i S i ~ H X^ReSi + Im5i ) (22) 



We find the same results as above, provided we identify the magnetization M with (— l/2)dF / dH . Other couplings 
to the external field or interpretations of the magnetization lead to quantitative differences. 

So far all models discussed have Gaussian properties since spins either are coupled in pairs or to an external field. 
However, in ref. [8] we have also allowed for additional random couplings between quartets of spins and possibly 
higher spin multiplets. It was found that in this situation replica symmetry is broken. Near the phase transition there 
appears a very close analogy with the SK model [9] for the mean field Ising spin glass. The replica symmetry breaking 
solution could be evaluated explicitly. The order parameter function has a unique shape at all T. In the classical 
description this model has again the unphysical low temperature effects of the spherical model. As, to our knowledge, 
it is the first spin glass model with infinite order replica symmetry breaking for which the whole low temperature 
phase has been solved exactly, it is an important question whether quantization using eq. (21) then leads a physically 
acceptable behavior in at low temperatures. It could indeed be shown that the entropy vanishes at T = 0. In fact, 
both entropy and specific heat again vanish as T 3 / 2 at low T. 

In conclusion, we have discussed a quantum description of the spherical model. For real spins we quantize Z 2 , for 
complex spins Z itself. The low temperature behavior is physical. Free spins in a field behave at low T indeed as 
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quantized objects. The same occurs in case of mean field couplings, since they lead to an effective field. For short 
range couplings the excitations are spin waves. In systems with uniaxial couplings the order will be suppressed in 
a large enough transversal field. The behavior near zero temperature quantum phase transitions is analyzed. Also 
uniaxial (Ising-like) and easy plane (XY-spins) systems can be described in appropriate limits. 

There occur two system parameters that can be adjusted to describe the correct Curie law at large temperatures 
and the correct zero point magnetization of real quantum spins. It is hoped that the spherical model can now be used 
in practical situations. 

ACKNOWLEDGMENTS 

This work is dedicated to the memory of Marc Kac. It is a pleasure to thank Bernard Nienhuis and Mark van 
Rossum for discussion. This work was made possible by the Royal Dutch Academy of Arts and Sciences (KNAW). 



[1] M. Kac, Phys. Today 10 (1964) 40 

[2] T.H. Berlin and M. Kac, Phys. Rev. 86 (1952) 821 

[3] G.S. Joyce, in Phase Transitions and Critical Phenomena 2, C. Domb and M.S. Green, eds. (Academic, London, 1972) 
[4] Y. Tu and P.B. Weichman, Phys. Rev. Lett. 73 (1994) 6 

[5] J.W. Negele and H. Orland Quantum Many- Particle Systems (Addison- Wesley, New York, 1988) 
[6] M. Kosterlitz, D.J. Thouless, and R.C. Jones, Phys. Rev. Lett. 36 (1976) 1217 
[7] Th.M. Nieuwenhuizen, Phys. Rev. B 31 (1985) 7487 

[8] Th. M. Nieuwenhuizen, 'Exactly solvable model of a quantum spin glass', preprint; cond-mat/9408056 
[9] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35 (1975) 1972 




0.0 1.0 2.0 



T/H 

FIG. 1. Specific heat of a free spin, a) quantum spherical model; b) classical spherical model; c) quantum 5 = 1 spin. 
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FIG. 2. Specific heat of S = 1 spherical spins on a simple cubic lattice: classical versus quantum situation. 
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AE/H 


1/2 


1/2 


0.23797 


0.96392 


0.95190 


1 


1 


0.41881 


2.25643 


0.83762 


3/2 


3/2 


0.59213 


4.02639 


0.78950 


2 


2 


0.76269 


6.28807 


0.76269 


5/2 


5/2 


0.93192 


9.04576 


0.74553 


3 


3 


1.10041 


12.30121 


0.73360 


7/2 


7/2 


1.26842 


16.05528 


0.72481 


4 


4 


1.43614 


20.30843 


0.71807 


9/2 


9/2 


1.60365 


25.06096 


0.71273 


5 


5 


1.77100 


30.31300 


0.70840 



Table caption: Parameters of the spherical model. 
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